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The superfluid-insulator transition of bosons is strongly modified by the presence of Fermions. 

Through an imaginary-time path integral approach, we derive the self-consistent mean-field transi- 
tion line, and account for both the static and dynamic screening effects of the Fermions. We find 
that an effect akin to the fcrmionic orthogonality catastrophe, arising from the fermionic screen- 
ing fiuctuations, suppresses superfluidity. We analyze this effect for various mixture parameters 
and temperatures, and consider possible signatures of the orthogonality catastrophe effect in other 
measurables of the mixture. 



I. INTRODUCTION 

The superfluid to insulator (SI) transition of bosons 
provides a conceptual framework for understanding 
quantum phase transitions in many physical systems, in- 
cluding superconductor to insulator transition in films 
wires®^^, Joscphson junction arrays^'^^^^, quantum Hall 
plateau transitions^^, and magnetic ordering^"*. The- 
oretical work on this subject elucidated many dra- 
matic manifestations of the collective quantum behav- 
ior in both equilibrium properties and out of equilibrium 
dynamics^^"^^. In many cases, however, we need to un- 
derstand the SI transition not in its pristine form, but in 
the presence of other degrees of freedom. For example, 
in the context of the superconductor to insulator transi- 
tion in films and wires, there is often dissipation due to 
fermionic quasi-particlcs, which may dramatically change 
the nature of the transition^^'^^"^^. Remarkable progress 
achieved in recent experiments with ultra-cold atoms in 
optical lattices (see ref. 23 for a review) makes these 
systems particularly well suited for examining quantum 
collective phenomena, not only as exhibited directly in 
the superfiuid phase, but also through its interplay with 
other correlated systems under study. 

A class of systems that can provide a new insight on 
the role of dissipation and of a fermionic heat bath on 
the superfiuid-insulator transition are Bose-Fermi mix- 
tures of ultra-cold atoms in optical lattices. Earlier theo- 
retical work on these systems focused on novel phenom- 
ena within the superfiuid phase, where coupling between 
fermions and the Bogolubov mode of the bosonic super- 
fiuid is analogous to the electron-phonon coupling in solid 
state systems. Several interesting phenomena have been 
predicted, including fermionic pairing^^^^^, charge den- 
sity wave order^^^^^, and formation of bound Fermion- 
Boson molecules^^''^^. Yet when Bose-Fermi mixtures 
were realized in experiments^^"^*, the most apparent ex- 
perimental feature was the dramatic loss of bosonic co- 
herence in the time of fiight experiments even for modest 
densities of fermions. This suggested an interesting pos- 
sibility that adding fermions can stabilize the Mott states 
of bosons in optical lattices. Theoretical work addressing 
these experiments, however, suggested that in the case of 



a homogeneous Bose-Fermi mixture at constant and low 
temperatures, the dominant effect of fermions should be 
screening of the boson-boson interaction, which favors 
the superfiuid state^^'^^'*°. Hence, the loss of coherence 
observed in experiments was attributed to effects of den- 
sity redistribution in the parabolic trap or reduced cool- 
ing of the bosons when fermions were added into the 
mixture. 

In this paper we argue that adding fermions into a 
bosonic system can actually stabilize bosonic Mott states 

even for homogeneous systems. While all previous the- 
oretical analysis represented the effect of fermions on 
bosons as an instantaneous screening, in this paper we 
take into account retardation effects, which arise from 
the presence of very low energy excitations in a Fermi 
sea. We show that such retardation gives rise to an 
eff(x-t which is analogous to the so-called orthogonality 
catastrophe, which is a well known cause for X-ray edge 
singularities and emission suppression^^'^^ in solid state 
systems. 

Our paper provides an alternative theoretical approach 
to the analysis of the Bose-Fermi mixtures. Rather than 
doing perturbation theory from the superfiuid state, we 
consider the Mott insulating state of bosons as our start- 
ing point. This is a convenient point for developing a 
perturbation theory, since deep in the Mott state the 
bosonic density is uniform and rigid and is accompanied 
by a simple Fermi sea of fermions. In general, the SI tran- 
sition can be understood as Bose condensation of parti- 
cle and hole-like excitations^'^ on top of a Mott state. In 
the absence of fermions, this condensation requires that 
the energy cost of creating particle and hole like exci- 
tations, i.e., the Hubbard U, is compensated by the ki- 
netic energy of these excitations, which is proportional 
to both the filling factor and the single particle tunnel- 
ing. Adding fermions to the system reduces the energy 
cost of creating either a particle or a hole excitation due 
to screening^^'^^, but it also reduces the effective tunnel- 
ing of bosons. The latter effect can be understood from 
the following simple argument. Consider a particle (or 
a hole) excitation on top of a Mott state of bosons. For 
fermions, this extra particle appears as an impurity on 
top of a uniform potential. When the bosonic particle 
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moves to the neighboring site, the "impurity potential" 
for all fermions changes. For individual fermionic states, 
the change of the single particle wave function may be 
small. But the effective tunneling of the bosonic particle 
is proportional to the change of the entire many-body 
fermionic wave function, and therefore we will need to 
take a product of all single particle factors. Even when 
each of the factors is close to one, the product of many 
can be much smaller than one. This is the celebrated 
"orthogonality catastrophe" argument of Anderson"'^. It 
can also be thought of as a polaronic effect in which tun- 
neling of bosons is strongly reduced due to "dressing" by 
the fermionic screening cloud. We see that both the inter- 
action and the tunneling are reduced by adding fermions. 
It is then a very non-trivial question to determine which 
effect dominates, and whether it is the superfluid or the 
insulating state that is favored by adding fermions into 
the system. Indeed, the main focus of our work is to un- 
derstand how the Fermi-Bose system pits the Bosonic su- 
perfluidity against the trademark dynamical effect of free 
Fermions. A related work, which addresses Fermionic dy- 
namical effects on the nature of the superfluid-insulator 
transition, is Ref. 44. 

In this letter we derive the SF-Mott insulator criti- 
cal line by constructing a mean-field theory that con- 
tains both the static screening and dynamical orthogo- 
nality catastrophe of the Fermions. For this purpose we 
resort to a new path-integral formulation of the mean- 
field Weiss theory for the SF-insulator transition^''. After 
demonstrating our approach by deriving the mean-field 
transition line for a pure bosonic system, we derive the 
path-integral approach to the Fermi-Bose system, and 
analyze the results in various limits. 

Our analysis will rely on several simplifying assump- 
tions. We consider only homogeneous systems, which 
is not the case for realistic systems in parabolic con- 
fining potentials. We do not allow formation of bound 
states between particles, which Umits us to small values 
of the Bose-Fermi interaction strength. The latter as- 
sumption becomes particularly restricting in one dimen- 
sional systems^^'^^"*'^, where even small interactions are 
effective in creating bound states. We do not take into 
account effects of non-equilibrium dynamics, which are 
important for understanding behavior of real experimen- 
tal systems whose parameters are being changed. And 
finally, we assume that there are only two fundamental 
states for bosons in the presence of fermions: the super- 
fiuid, and the Mott insulator. When our analysis shows 
proliferation of particle and hole like excitations inside 
a Mott state, we interpret this as the appearance of the 
superfluid state. We do not consider the possibility of ex- 
otic new phases such as the compressible state suggested 
recently by Mering and Fleischhauer'*''' . While these limi- 
tations make it difficult to make direct comparison of our 
findings to the results of recent experiments^^"^^, we be- 
lieve that our work provides a new conceptual framework 
which can be used to address real experimental systems. 



A. Microscopic Model 

The Hamiltonian for the Bose-Fermi system we analyze 
is given by 

H = Hb + Ti-F + "Hint 
= Elif^s^i - l^ni-lY, J cos - <i)j)] 

i {ij) 

Tip = -J2 Jpclcj -Y,iiFc\ci, Hint = YlUpBriiclci. 

{ij) i i 

(1) 

Hb describes the bosonic gas using the phase and num- 
ber operators in each well: hi, (pi. J is the strength of the 
Josephson nearest-neighbor coupling (note that J « nt, 
where n is the filling factor and t is the hopping amplitude 
for individual bosons) , Ub and fi are the charging energy 
and chemical potential, respectively. Hp describes the 
Fermions, with hopping Jp and chemical potential Hf- 
The two gasses have the on-site interaction UpB- For 
simplicity we use the rotor representation of the Bose- 
Hubbard model, but our results are easily generalized 
to its low-filling limit. The pure Bose gas forms a su- 
perfluid when J/A ~ 1, where A ~ J7 is the charging 
gap^'''^^. The fermions encourage superfluidity, on the 
one hand, by partially screening charging interactions 
and reducing the local charge gap^^'^°. But at the same 
time the fermions' rearrangement motion in response to 
boson hopping is slow and costly in terms of the action it 
requires. This motion results in an orthogonality catas- 
trophe that suppresses superfluidity. 

Our derivation of the phase diagram is based on the 
mean-field approach, which in the case of purely bosonic 
systems is equivalent to the analysis in Refs. 16,17,48, 
but which can be generalized to study Bose-Fermi mix- 
tures. The idea is to use the Weiss approach of reducing 
the many-site problem in the Hamiltonian (1) to a single 
site problem by assuming the existence of the expectation 
value for the phase coherence of bosons: 

(e^^O=r-. (2) 

In the local problem, one can calculate a self-consistent 

equation for r that will produce the transition point. 
This procedure can not be simply followed once the 
fermions are thrown into the mix, since even with Eq. 
(2), the Hamiltonian H is non-local; this problem is ad- 
dressed by using the imaginary-time path integral formu- 
lation. 



B. Overview 

This paper is organized as follows. In Sec. II we derive 

the path-integral formulation for the mean-field phase 
boundary of a pure Bose g function of the param- 

eters in its Hamiltonian and temperature. In Sec. HI 
we build on this formalism to account for the weakly in- 
teracting Bose-Fermi mixture. We find a new condition 
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for the supcrfluid insulator transition in terms of the Bo- 
son parameters, as well as the interaction strength, Ufb, 
and the Fermion's density of states, p. Our main result 
is presented in Sec. IV in Eq. (30). The mean-field con- 
dition is plotted for the cases of fast and slow Fermions, 
for zero temperature, as well as at a finite temperature. 
We conclude the paper with a summary and discussion 
in Sec. V. 

Our main findings are that even a moderately weak in- 
teraction with slow Fermions inhibits superfiuidity in the 
Bosons. The dynamical response of the slow Fermions 
produces a large cost in terms of the action for bosonic 
number fiuctuations. This effect of the orthogonality 
catastrophe of a Fermionic screening gas is most appar- 
ent where the on-site charging gap of the Bosons is small 
(A <C U). In Sec. IV we also derive approximate simple 
expression for the phase boundaries for this case at zero 
and low temperatures, Eqs. (31) and (32). Our analysis 
shows that the phase boundary becomes non-analytical, 
and superfluidity is dramatically suppressed. 



equation, Eq. (2), becomes: 

r = ^ J D[<j>jiT)] E cos<^,(0)exp(-5,) 

^r^jD[<f>{T)] E /dricos</.(0)cos(.^(Ti)) 

Mr)} 



(6) 



exp 



J dr —in^+^Un^ — fin 
'- 



where cos ^(0) was expanded in r to its lowest power. 

Our goal is to simplify condition (6); for this purpose, 
we integrate over the phase variable (p. Let us concentrate 
first on the partition function, Z, in the denominator 
of Eq. (6). In the limit of r — > 0, </> only appears in 
the Berry-phase term, which using integration by parts 
becomes: 

/3 

i / dTn4) = in(0) ((/)(/?) - (j){0)) - i / dr^n. (7) 



II. PURE BOSE-GAS PHASE-DIAGRAM 
USING THE PATH-INTEGRAL APPROACH 

We begin our analysis with the pure bosonic gas. We 
will use this case, where no Fermions are present, to de- 
rive and demonstrate our path-integral approach to the 
mean-field superfluid-insulator transition. We will first 
use the mean- field ansatz, Eq. (2), to reduce the parti- 
tion function to a path-integral over a single-site action. 
Analyzing the single-site action will reveal the mean-field 
condition for superfluidity. 

The flrst step is to transform the Hamiltonian Hb of 
Eq. (1) into a single-site Hamiltonian. Using Eq. (2), we 
can write: 



Because n(0) is an integer and (/)(r) is periodic on the 
segment [0, the flrst term is always a multiple of 27ri, 
and can be omitted. Furthermore, without an r term, the 
angle variables in each time slice, ^(t), become Lagrange 
multipliers which enforce number conservation in the site: 



j D[(l){T)]exp{-i j dr^Ti) = Jj27r^^(^),o = Jj27r^„(^),„ 

(8) 

Thus n(r) = 0, and we can write the single-site partition 
function as: 



(0)- 



(9) 



1, 



Kb ^ HBj = ^Ubu] - fih, - z-J {r*e"l'' + re"*"^^) 



(3) 

where z is the coordination of the lattice. The action for 
site j therefore becomes: 



Sj = I dr 



Thus the partition function for a single site is: 

— / dr^i{pn— zJr cos 4'+^ Un'^ — tin^ 



Z = J D[</.(r)] 



e 



(4) 



(5) 



Mr)} 



where we assumed that r is real, and dropped the index 
j- 

The self-consistent condition for superfluidity equates 

the degree of phase ordering on site j with r, which was 
substituted for the neighbors of site j. This mean held 



Eq. (6)'s numerator is more involved. The phase 0(t) 
now also appears through the cos(0(ti)) cos(0(O)) term. 
A e^'^^'^i) term is indeed a creation operator, therefore we 
expect that the cosine factors in the path-integral will 
change the number of particles n(r) at ri and at r = 0. 
Let us demonstrate this by concentrating on the term 
cos?!>(0) = i(e''^(°' + e-'-^o) and integrating over ^(0): 

TT ('^An(0),l + '^An(0),-l) ^ 

where An(T) ~ h(T)dT ~ lim£^o'^(''' + ^ '^(''" ^ f ) ■ The 
same expression results from the integration on the ti 
time slice. Thus the integration of the 4>{t) variables in 
the numerator of Eq. (6) still gives n(r) = as long as 
r 0,Ti. The numerator of Eq. (6) can now also be 
reduced to a simple sum over n(T) , but with a jump in 
the boson- number at tq and ri: 



n(r) =n±e{T)0{Ti -t). 



(11) 
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Now that the integration over the (/)(r) variables is 
complete, we can write the mean- field condition for su- 
perfluidity as a single sum. Since the -I- and - choices 
in Eq. (11) give rise to the same contribution in the 
mean-field condition, we can choose the plus, n(r) = 
n + 6{t)9{t\ — r), and write: 



/3 
/ 

-oo 







n=—oo 



5 exp 



— / drHcin) 



(12) 



exp[- / drHcin + 1)], 





with Hcin) = ^JJ-n? — fin. For a pure Bose gas, we obtain 
the well known Weiss mean-field rule for X-Y magnets: 



1 = 



zJ/2 



OO 



,-/3W,.(n) 



,-/3H,(n+l) 



oo ^ ^ 

— — oo 



Hc(n+1) -Hc(n) 



(13) 



III. EFFECTIVE BOSONIC ACTION FOR THE 
SF-INSULATOR TRANSITION OF THE 
FERMI-BOSE MIXTURE 



By finding the minimum with respect to the Fermion 
density, uf, we find: 



nF = n%- UpBpnB, 
and the total charging energy is: 



(15) 



E^ = Eo+-{U- UIbp) nl-{iiB - UfbuV) ub- (16) 

Therefore the charging parameters of the Bose gas are 
renormalized by the presence of the Fermions to: 

U=U-UIbP ilB = IJiB-UFBn%. (17) 

This charging energy renormalization makes the Mott 
lobes shrink in the p-J parameter space by the ratio U /U: 
adding the Fermions mitigates any static charging effects, 
since the mobile fermions c;an screen any local charge 
even when the bosons are localized. 

An important note is that in order for the hydrody- 
namic approach to be correct, the electronic screening 
should not exceed one particle. This restricts the pertur- 
bative regime to: 



The addition of Fermions to the Bosonic gas affects 
the bosons in two distinct ways. The first is static: the 
Fermions shift the chemical potential and the interac- 
tion parameters of the Bosons^^'^^. But in the superfiuid 
phase. Boson number fluctuations become dominant, and 
the screening problem becomes a dynamical one. The 
Fermi screening cloud requires a flnite time to form, and, 
in addition, it costs an prohibitively large action in some 
cases. While the former static screening effect enhances 
superfluidity, the latter dynamical effect suppresses it. 
The advantage of the imaginary-time path integral for- 
malism, which was developed in the previous section, is 
that it deals with both effects on the same footing, and 
allows the inclusion of the fermionic collective dynamical 
response in a one-site bosonic action. 



A. Static screening effects 

Let us now consider the Fermi-Bose mixture of Eq. 
(1). The most straightforward effect of the Fermions is to 
shift the chemical potential and interaction parameters. 
We will first calculate this effect using a hydrodynamic 
approach^^'^^. By denoting the DOS of the fermions at 
the Fermi-surface as p, and neglecting its derivative, we 
can write a charging-energy equation for the mixture per 
site: 



with 



E\^\nF) - i^ftIf +^Un% - pbtib + UFBnBnp 

(14) 



dn 



Ep 



dni 



UfbP < 1- 



(18) 



In addition, for the Fermi-Bose mixture to be stable, we 
must have U > 0, and thus also: 



U'fbP < U. 



(19) 



which in the regime of interest is a less restrictive con- 
dition than Eq. (18). Next, we consider the Fermionic 
dynamical response. 



B. The Fermion's dynamical response 

The superfiuid bays between the Mott lobes in the tra- 
ditional fi — J phase diagram are affected strongly by a 
more subtle and intriguing effect: dynamical screening 
motion of the Fermions. The analysis of this effect makes 
the path integral necessary. We construct the path inte- 
gral starting with the action: 



Sfb = J dr l-i E(CiCi -|- UBi'pi) + Hf + 'Hb+ 
E Ufbu^ubi + UfbUbi {c\ci - n^^ j 



(20) 



with -nPp = (cjci)) and where c and should be con- 
strued as Grassman variables. The first term in the sec- 
ond line produces the shift in the chemical potential, as 
in Eq. (17), but U is not yet renormalized. The U 
renormalization is a second order effect, which we analyze 
by producing a perturbation series in UpB^np i, where 



Arif 



An effective action for the Bosons 
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is then obtained by integrating over the Fermionic vari- 
ables: 

e °B = e ' 

- fdrl -iJ2{cjCj)+HF + UFBnBj^nFj ) 
Tie V J / Rj 



jD[c\D[c^]e 



Zpe 



ghUpB J dn J dT2nj{Ti){AnF j{ri)AnF j{T2))n.i{T2) 

(21) 

where Tis is the pure bosonic Hamiltonian with the 
renormalized charging energy and chemical potential, Eq. 
(17). 

From Eq. (21) wc see that the integral over the 
Fermionic degrees of freedom gives rise to a new Boson 
interaction term. It is given as a polarizability bubble for 
the fermions: 



fei, k2 



5 



(22) 



with 



Hf being the fermionic kinetic energy 



relative to the Fermi surface. After some manipulations 
(see App. A) we obtain: 



screening cloud forms only after the time p. The instan- 
taneous screening is thus modified by the action (when 
Tp«l): 



-(ti — y9 In [cosh - 



cosh - 



- cosh - 



cosh 



2 



'-])■ 



(26) 



While Scharging assumcs that the polarizability, Eq. (23) 
has the term ap for all frequencies, the correction term 
takes into account the cutoff in the static screening 
term. Instead of the fermionic screening in the wake of 
a change of at r = being An_F ~ ^('''), we have 
Arij^ ^ arctanh(T//7). Considering in addition the peri- 
odic nature of imaginary time, we obtain Eq. (26). 

Finally, SqCi contains the contribution due to the or- 
thogonality catastrophe: 



ttU^bp'' , sin(7rrTi) 



7ln — 



sin(7rT/A) ' 
(27) 

In Eq. (27) we defined the dissipation parameter: 



(28) 



Note that we are restricted to 7 < 1 in the perturbative 
regime, Eq. (18). 



{AriFiAnF iiO))u 



ap-Tr\u}\p \u)\< p 1 

c/a;2 \uj\ > 1/p, ^^•^> 



and as we discuss below, the perturbative analysis is valid 
when pUfb < 1- The first low- frequency term in Eq. 
(23) yields the static screening, Eq. (17), i.e., a = Up^p. 
The second term describes the dynamical response, pro- 
ducing the action term: 



(24) 



This term yields logarithmic contributions to the ac- 
tion, whose effects are familiar from electronic systems 

as orthogonality catastrophe in metal X-ray absorption 
spectrum, the Kondo effect^", and Caldeira-Leggett 
dissipation^^'^^. Here its effect is to suppress superfluid- 
ity, since it couples to the number fluctuations. At angu- 
lar frequencies greater than 1/p, the interaction term de- 
cays quickly (c is a positive constant), and hence A = 1/p 
serves as a UV cutoff. It also implies that the screening 
and logarithmic contributions to the action can only ap- 
pear with a time-lag Td ^ p. 

The next step is to calculate the action <S'[(n(r)] with 
n(r) =n + 9{t)0{ti - t) as in Eq. (12). We have: 



S[n{i 



^charging 



^oc 



(25) 



First, Scharging = Jq i^Uni^T)'^ - pn{T)) . But the cutoff 
A « 1/p of the fermions' polarizability implies that the 



IV. MEAN-FIELD PHASE DIAGRAM AND 
THE ORTHOGONALITY CATASTROPHE 

The mean-field transition line is obtained, as in Sec. 
II, from Eq. (12), which here takes the form: 



zj_ 
2Z 



= -S[n(r)] 



(29) 



n=— 00 



Substituting S[n{T)] = Scharging + from Eqs. 

(26) and (27), we obtain the mean-field condition for the 
transition line: 



Z 



exp[-r(We(n + 1) - H.{n)) - S\t)] • (i^g^)' , 

(30) 

with Hc{n) = ^Un^ — fin, and S^{t) given in Eq. (26). 
This is our main result. 

Eq. (30) allows us to calculate the mean-field SF- 
insulator phase boundary for weakly interacting mixtures 
for a range of Temperatures and Fermi DOS. To illus- 
trate Eq. (30) predictions for the transition line. Figs. 
1, and 2 show the boundaries for Bosons and Fermions 
interacting with Ufb = O.IU and Ufb = 0.25U, respec- 
tively, for a range of fermion velocities, or DOS p. In Fig. 
3 we plot the effect of slow fermions, with p = 3.5/U, 
on the Bosonic SF transition for a range of interactions 
Ufb- The Superfluid-insulator transition boundary at 
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FIG. 1; Zero-temperature mean-field phase diagram of the 
Fermi-Bose mixture with Ufb = 0.1J7. (a) Renormal- 
ized chemical potential, jl/U vs. bare J/U phase bound- 
ary. From bottom-left to top-right, the fermions DOS is 
p ■ U — 3.5, 3, 2.5, . . . , 0.5. A dashed red line marks the un- 
coupled Bose gas, Ufb ~ 0, but can barely be distinguished 
from the p = 0.5 line, (b) charge-offset n = p,/U vs. bare 
J/U phase boundary. From left to right at n = 0, the fermions 
DOS is p- (7 = 3.5, 3, 2.5,..., 0.5. 



finite temperature for UpB — 0.25U is shown in Fig. 4 
for a range of temperatures for fast and slow fermions. 
In all figures we assume z = 6. 

One easily drawn qualitative conclusion is that slow 
electrons mostly inhibit superfluidity, which is the mark 
of the orthogonality catastrophe. The most dramatic 
suppression effect occurs near the degeneracy points, 
where n — jl/U = m + 1/2. Let us obtain closed-form 
expressions for the SF-INS boundary there. A helpful 
observation is that if the charging gap nearly vanishes, it 
suffices to consider the lowest nearly-degenerate charge 
states in Eq. (30). 

When the degeneracy is exact, e.g., at fT = 1/2, we 
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FIG. 2: Zero-temperature mean-field phase diagram of the 
Fermi-Bose mixture with Ufb = 0.25f/. (a) Renormal- 
ized chemical potential, p,/U vs. bare J/U phase bound- 
ary. From bottom-left to top-right, the fermions DOS is 
p-U = 3.5, 3, 2.5, . . . , 0.5. The dashed red line is the un- 
coupled Bose gas, Ufb = O (b) charge-offset n — jl/U vs. 
bare J/U phase boundary. From left to right at n = 0, the 
fermions DOS is p • [/ = 3.5, 3, 2.5, . . . , 0.5. (c) A focus on 
the area near degeneracy, n = 1/2, in (b). In this plot we can 
see the strong effect of the orthogonality catastrophe for slow 
electrons. The intercept of the boundary with the n axis be- 
comes singular and scales as J/U ~ |0.5 — n\^~'' , illustrating 
Eq. (32). 
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FIG. 3: Zero-temperature mean-field phase diagram of the 
Fermi-Bose mixture with Fermi dispersion p = 3.5/(7. (a) 
RenormaUzed chemical potential, jl/U vs. bare J/U phase 
boundary. From bottom-left to top-right, the fermions DOS 
is Ufb/U = 0,0.05, 0.1,..., 0.25. Ufb = is shown as a 
dashed red line, (b) Charge-offset n — jl/U vs. bare J/U 
phase boundary. From left to right at n = 0, the fermions 
DOS is Ufb/U = 0.25, 0.2, . . . , 0, with Ufb = shown as a 
red dashed line. 



obtain for the critical J vs. Temperature: 





-1 zj 

~ 4 



zJ ( 2TrT \1 
AttT V a / 



( sn,(.T/A) y 



Tr) J 

7) sin I 



(31) 



£7 1 
2 J 



where we neglected screening retardation, i.e., 5*^, alto- 
gether. This is valid for T^-t < B{m,n) is 
the beta-function. Fig. 4b shows Eq. (30) in this limit. 
Whereas for pure bosons the critical J is linearly pro- 
portional to T, the interaction with Fermions makes the 
critical J required for superfluidity increase dramatically, 
and obey Ja ~ T^^t. 

A similar analysis can be done at zero temperature 
slightly away from degeneracy atn — l/2 = e^l. In 




Ufb=0.25 
^002 004 006 OM 



T/U 



FIG. 4: Temperature dependence of the SF-Ins boundary 
for Ufb = 0.25!7. (a) The SF-Ins boundary for T = 
0,0.04, ... ,0.2 (from left to right at n = 1/2) with fast 
fermions, p = 0.5. The dashed red line is the uncoupled Bose 
gas boundary, Ufb = 0. (b) Same as (a) with slow Fermions, 
p = 3/(7. (c) The critical boson hopping J/U vs. T/U, for 
fermion DOS p ■ U = 0.5, 1, . . . , 3.5. The linear curve of the 
uncoupled Bose gas becomes a cusp as the bosons couple to 
slow fermions, as also calculated in Eq. (31). 
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this regime we obtain: 

^ ^ 

where ignoring screening retardation is vaUd if e^~^ <^ 
{K/UY~^ Here too, the critical hopping as a function 
of e is linear, Jc ^ |e| for pure Bosons, but increases 
dramatically to kl^ when the Bosons interact 

with the Fermions. This dependence is demonstrated in 
Fig. 2c. 



Experiments on Bose- Fermi systems show a strong sup- 
pression of superfluidity. Ref. 32 describes a system 

where J/U ~ 1/20, Jp/J ~ 5, and UpB 2C/, i.e., 

the Bose-Fermi system is strongly interacting. Thus our 
theory of orthogonality-catastrophe effects is not directly 
applic;ablc here. We note, however, that at large val- 
ues of UpBi bound composite fermions would form^°. 
These will have a weakened Ups and a strongly enhanced 
DOS, p. This makes it possible to observe orthogonality- 
catastrophe SF suppression even in this regime. 



V. DISCUSSION 

A. Regime of applicability 

Our theory of the SF-insulator transition applies to 
the weakly coupled Bose-Fermi mixtures, with Ufb < 
U, but a second parameter that is required to be small 
is UfbP < 1 (see Eq. 18). The latter is required for 
the perturbation theory of Eq. (21) to be justified. As 
explained above, this condition can be easily understood 
by noting from Eq. (15) that the response of the Fermi 
gas to the appearance of a Boson in a particular site 
is Anp = UfbP; which clearly must be lower than 1. 
Even more importantly, the perturbative analysis is valid 
so long that no localized states form in the fermionic 
spectrum when a site's potential changes by UpB' this, 
too, is true when pUpB < 1 at large dimensionality. 

The formation of a localized state at larger values of 
Up BP, and therefore where 7 > 1, is likely to suppress 
the orthogonality-catastrophe effects, perhaps in analogy 
to the behavior of a Kondo-impurity in a metal: when a 
Kondo impurity localizes an electronic state it becomes 
inert. Therefore the largest suppression of the SF-INS 
boundary is likely to occur when UpbP ~ 1, as 7 ~ 1^^. 
The regime 7 > 1 lies beyond the scope of this paper, 
but we intend to approach it in a later publication. Note 
that this regime can still occur when UpB ^ U . 

We note also that since our theory is concerned with 
weak Fermi-Bose coupling, it ignores Fermi-Bose bound 
pair formation, as well as p-wave superconducting corre- 
lations, which may be important only at parametrically 
low temperatures. 

B. Relation to experiment 

Our theory provides the mean-field phase diagram un- 
der the assumption of a grand-canonical ensemble with 
fixed chemical potentials. Experiments, on the other 
hand, are conducted in finite non-uniform traps, and 
therefore to compare their results with the thermody- 
namic phase diagram we provide, the chemical potential 
of the interacting Bose and Fermi gasses must be de- 
termined for particular trap geometries, using, e.g., the 
LDA approach, as in Refs. 24,25. 



C. Relation to dissipative phase transitions 

At weak UpB, we find that Fermions, through their 

orthogonality catastrophe, by and large inhibit superflu- 
idity, particularly when the fermions are slow. This effect 
is extremely reminiscent of clissipative superconducting- 
metal phase transitions in Josephson junctions. 

Typically, dissipative effects as in resistively shunted 
Josephson junctions (RSJJ), are thought to strengthen 
phase coherence^'*'^^. But in our case, since the Bose- 
Fermi mixtures couple through a capacitive interaction, 
as opposed to the phase-phase interaction in super- 
conducting systems expressed in a Caldeira-Leggett"''^''''^ 
term or its modular equivalents^, we encounter a sup- 
pression. 

Another important distinction is that the Caldeira- 

Leggett analysis of a single RSJJ's, and of Id superflu- 
ids, associates (Quasi long-range) phase ordering with the 
long-time behavior of (e*'^('^)e~"^*^°^). But in the mean- 
field theory of the SF-insulator transition, the onset of 
true long-range phase order we encounter is associated 
with the less restrictive time-integral of the aforemen- 
tioned correlation, as in Eq. (30)^^. The dynamics of 
the Fermionic screening gas modifies this integral only 
quantitatively, but it does not affect the nature of the 
transition. 

The less restrictive condition for ordering in the mean- 
field analysis reflects the assumed higher dimensional- 
ity of the systems we consider. Concomitantly, in low 
dimensional systems with short-range interactions only, 
the Mermin- Wagner theorem rules out the formation of 
long range order. For Josephson junctions, for instance, 
phase-slips are the domain-wall Uke defects which make 
the order parameter fluctuate. Similar defects are ab- 
sent from our analysis since their cost in terms of action 
is too prohibitive due to the assumed high connectivity 
of the system. Therefore we can make the mean-field 
assumption of a non-fluctuating order parameter. This 
assumption is fully justified above the lower critical di- 
mension (although our analysis will only be valid at and 
above the upper critical dimension). 
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D. Summary and future directions 

In this manuscript we concentrated on the effects of 
the orthogonality catastrophe on the superfluid-insulator 
transition line, and showed how slow fermions inhibit 
superfluidity through dissipation capacitatively coupled 
to number-fluctuations. The orthogonality catastrophe 
should also be evident in other measurements, which may 
give an independent estimate of the dissipation present. 
This might be most apparent in revival experiments, 
where the system is shifted from a superfluid phase into 
the insulating phase, and released after t^^,^^. We expect 
that the revival decay time will be smaller with increasing 
dissipation. We intend to address the dynamical effects 
in Fermi-Bose mixtures in a future work. 



Another interesting angle for future work is the ap- 
pearance of a super-solid at the special point of Fermionic 
half-filling^^, extending our formalism to account for this 
possibility could be done by considering the Fermionic 
density correlations near nesting vectors of the Fermi gas. 
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APPENDIX A: FERMIONIC RESPONSE FUNCTION 



For completeness, we provide here a simple derivation of the Fermionic response function, Eq. (23), for fermions 
in one-dimension. Once our result is put in terms of the fermionic density-of-states at the Fermi- surface, it applies 
in any dimension, since the existence of a {d — l)-dimensional Fermi-surface renders the dispersion for low-energy 
excitations essentially one-dimensional. 

Let us assume, for simplicity, that the fermionic Hamiltonian is: 

2kF 

n = hVF - ^F)(4cfe + clfcC_fe) (Al) 

fe=0 

The density of states per site for this Hamiltonian is: 



For the Hamiltonian (Al), Eq. (22) reads: 



1 ''^ -1 

This formula sums over the contributions of particle-hole excitations of four kinds: both particle and hole are right 
movers {ki,k2 > 0), both particle and hole are left movers (fci,fc2 < 0), and two mixed cases. In order to avoid the 
absolute value, we concentrate on the first case, and multiply by four: 

C =4T Y r ^ r ^ (A4) 

where we also shifted w' by w/2. We now separate from this sum the contributions from large > fivrkp. These 

high energy modes contribute an w-independent term to the static screening. Corrections to this constant are easily 
seen to be quadratic in uj (e.g., set fci, ^2 0). The low-w' terms, however, will give rise to a contribution, which 
we are after. The k-integrals can be easily done in the limit \co' ± a;/2| <C fivrkp, and we obtain: 

n — — AT \^ 1 1r,c- f ^fFfcF-»(t^'-t^/2) \ 1 / hvpkF-i{uj'+ul/2) \ 

a,'=27r(n-M/2)-Fa,/2 ' ^ ^ /V V ^ 'U ^^^^ 

^c-AT E (2^^ [-iT^ sgn(u;' + a;/2)] [-Itt sgn(a;' - . 

bj'=2-K(n+l/2)+uj/2 



The w dependence arises from the region where the two sign functions give opposite results: —\uj\/2 < w' < |w|/2. 
Thus: 
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where the last factor of 2 is since d contains the contributions for — |a;|/2 < w' < |w|/2 assuming the same sign as for 
the rest of the frequency range. The final answer is thus; 



as reported in Eq. (23). 



1. Bosonization approach to the polarization calculation 



One-dimensional fermionic systems are most effectively described in terms of a bosonized action. Let us re-derive 
Eq. (23) using this simpler approach. We define the two fields d and Using here the convention: 

-Ve = pl^ PR -V^ = PR- PL, (A8) 

TT TT 

where pr^l are the right-moving and left-moving densities respectively, the Hamiltonian of 1-d Fermions is: 

n='^jdx{ve' + vci>'), (A9) 

and t;ir is the Fermi velocity. 

The Hamiltonian (A9) can be turned into an imaginary time Lagrangian: 

C='^ j dx {V9^ + 9^) . (AlO) 

The density-density correlation we would like to calculate is now given as a path-integral over the field: 

= (An^-,AnF,(0))„ = 2 • \{veve)^ (All) 

Here we need to pause and explain the extra factor of 2: the expectation value in the brackets only takes into account 
particle-hole excitations that are contained within the same branch of the fermionic spectrum, right moving or left 
moving. We must also include, however, excitations with the partic;le part being a right mover and the hole being 
a left mover, and vise versa. These give exactly the same contribution (as is also seen in the first approach), and 
therefore it is sufficient to simply multiply the expectation value by 2. 
With that in mind, we proceed to write: 

_ 2^ r dk k^ 



kp- 

f dk 

— kp 



(A12) 



2^ J 



dk 

fe2 + (w/vF)2 



= c'- = c'-7r/,2|^|, 

where in the last step we simply calculated the residue of the k integral. 
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